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Consider an x n random matrix Yn = (Y^^) where the entries 
are given by Y^" = "'^^ X"j , the X," being independent and iden- 
tically distributed, centered with unit variance and satisfying some 
mild moment assumption. Consider now a deterministic N xn matrix 
An whose columns and rows are uniformly bounded in the Euclidean 
norm. Let E„ = y„ + An. We prove in this article that there exists a 
deterministic N x N matrix- valued function Tn{z) analytic in C — R"*" 
such that, almost surely, 

lim ( ^Tr&ce{T.nT.l - zInV^ - ^Tr&ceTn{z)) =0. 

Otherwise stated, there exists a deterministic equivalent to the empir- 
ical Stieltjes transform of the distribution of the eigenvalues of E„E^. 
For each n, the entries of matrix Tn{z) are defined as the unique solu- 
tions of a certain system of nonlinear functional equations. It is also 
proved that Trace Tn{z) is the Stieltjes transform of a probability 
measure 7r„(dA), and that for every bounded continuous function /, 
the following convergence holds almost surely 

T7V/(Afe)-/ /(AK4dA)-^0, 

fc=l 

where the (Afe)i<fe<jv are the eigenvalues of E„S^. This work is 
motivated by the context of performance evaluation of multiple in- 
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puts/multiple output (MIMO) wireless digital communication chan- 
nels. As an application, we derive a deterministic equivalent to the 
mutual information: 

C4a') = lElogdet(/^ + , 

where is a known parameter. 



1. Introduction. 

The model. Consider an x n random matrix Yn where the entries are 
given by 

(1.1) y.';=^^?r 



n 

where (crjj(n), 1 < i < A^, 1 < j < n) is a bounded sequence of real num- 
bers (i.e., sup„max(jj) |crjj(n)| = Cmax < +oo) cahed a variance profile and 
the real random variables Xfj are centered, independent and identically dis- 
tributed (i.i.d.) with finite A + e moment. Consider a real deterministic N xn 
matrix An = (A'^'j) whose columns (a^)i<fc<„ and rows {sL^)i<i<N satisfy 

(1.2) supmax(||a^||,||a^||) <+oo, 

n>l 

where || • || stands for the Euclidean norm. Denote by T,n = Yn + An- This 
model has two interesting features: the random variables are independent 
but not i.i.d. since the variance may vary and An, the centering perturbation 
of Yn, can have a very general form. Denote by Tr(C) the trace of matrix 
C. The purpose of this article is to study the behavior of 

^Tr(S„S^-z/^)-\ ZGC-M+, 

that is, the Stieltjes transform of the empirical eigenvalue distribution of 
when n — > -|-oo and N +oo in such a way that ^ ^ c, < c < +00. 
Under Condition (1.2), the convergence of the empirical distribution of the 
eigenvalues of may fail to happen (see, e.g.. Section 3.1). We shall 

prove that there exists a deterministic matrix- valued function Tn{z), whose 
entries are defined by a system of equations, such that 

4Tr(S„S^-z/^)-^-4Trr„(z) 0, zGC-M+. 

N N n-*oo,N/n-*c 

It is also proved that -^TrT„(z) is the Stieltjes transform of a probability 
measure 7r.„((iA) and that the following convergence holds true: 



(1-3) j;^i:nx>^)-l mMdX)^-^^o, 
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where / is a continuous and bounded function and where the {Xk)i<k<N 
are the eigenvalues of Sn^^n- The advantage of considering j^TtTniz) as 
a deterministic approximation instead of E-^Tr(S„,S^ — zIn)~^ (which is 
deterministic as weh) hes in the fact that Tn{z) is in general far easier to com- 
pute than E-^ Tr(S„S^ — z/at)^^ whose computation relies on Monte Carlo 
simulations. These Monte Carlo simulations become increasingly heavy as 
the size of matrix S„, increases. 

Motivations. This problem is mainly motivated by the context of perfor- 
mance analysis of multiple input /multiple output (MIMO) digital communi- 
cation systems. The performance of these systems depends on the so-called 
channel matrix whose entries {Hfj, l<i<:N,l<j<n) represent the 
gains between transmit antenna j and receive antenna i. Matrix Hn is of- 
ten modeled as a realization of a random matrix. In certain contexts, the 
Gram matrix HnH* is unitarily equivalent to a matrix {A^ + 1^,)(^„, + Yn)* 
where Yn is a random matrix given by (1.1) and An is a possibly full rank 
deterministic matrix satisfying (1.2). A fundamental result in information 
theory (see [26]) states that the mutual information Cn, which is a basic 
performance index of the channel, can be expressed in terms of the singular 
values of the channel matrix Hn ■ 



where cr^ represents the variance of an additive noise corrupting the received 
signals. We shall show in Section 4 how to approximate such a mutual in- 
formation with the help of the deterministic equivalent Tn{z). Related work 
can be found in [7, 22, 26]. 

About the literature. If Z„ is a zero- mean N x n random matrix, the 
asymptotics of the spectrum of the N x N Gram random matrices ZnZn 
have been widely studied: Marcenko and Pastur [19], Yin [28], Silverstein et 
al. [5, 24, 25] for i.i.d. entries, Girko [9], Khorunzhy, Khorunzhenko and Pas- 
tur [17], Boutet de Monvel, Khorunzhy and Vasilchuk [4] (see also Shlyakht- 
enko [23]) for non-i.i.d. entries. In the centered case, it turns out that very 
often the empirical distribution of the eigenvalues converges toward a limit- 
ing distribution. 

The case where matrix Z„ has nonzero mean has been comparatively less 
studied. Among the related works, we first mention [5] which studies the 
eigenvalue asymptotics of the matrix {Rn + Yn){Rn + ^n)^ in the case where 
the matrices Yn and Rn are independent random matrices, Yn has i.i.d. 
entries and the empirical distribution of RnRn converges to a nonrandom 
distribution. It is shown there that the eigenvalue distribution converges 
almost surely toward a deterministic distribution whose Stieltjes transform 
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is uniquely defined by a certain functional equation. The case {Yn + A„) {Yn + 
An)'^ where Yn is given by (1-1) and A„ is deterministic pseudo-diagonal 
has been studied in [13] where it is shown that under suitable assumptions 
the eigenvalue distribution converges. In the general case S„ = Yn + An, the 
convergence of the empirical distribution of the eigenvalues of may fail 

to happen even if the variance profile exists in some limit and the spectral 
distribution of AnAn converges (see, e.g., Section 3.1). 

Girko proposed in [10], Chapter 7, to study a deterministic equivalent of 
some functional of the eigenvalues of S„E^ in the case where the following 
condition holds for An'. 

n N 

(1.4) supmaxV] I A" I < +00 and supmaxV] [^4" •[ < +oo. 
" ' 3=1 ' ^ ^ i=l 

Girko showed that the entries of the resolvent (S„S^ — zl)~^ have the same 
asymptotic behavior as the entries of a certain deterministic holomorphic 
N X N matrix- valued function Tn{z) characterized by a nonlinear system 
of (n + N) coupled functional equations. Condition (1.4) is however rather 
restrictive. In particular, it does not hold in the context of wireless MIMO 
system in which (1.2) is actually much more relevant. In this paper, we thus 
extend some of the results of Girko [10] to the case where An satisfies (1.2). 
For this, we do not follow the approach of Girko [10] based on the use of 
Cramer's rule but rather take the approach of Dozier and Silverstein [5] as 
a starting point. This approach not only allows to extend the result of [10] 
to deterministic matrices satisfying (1.2), but also provides a simpler proof. 
It also enables us to prove that the deterministic equivalent -^TrT'„(z) of 
the Stieltjes transform ■^Tr(S„S^ — zl)~^ is itself the Stieltjes transform 
of a probability measure, which is a result of interest from a practical point 
of view [see (1.3), e.g.]. 

Outline of the paper. The main results of the paper are Theorem 2.4 
[existence of the deterministic matrix-valued function T„(z)], Theorem 2.5 
(deterministic approximation of the Stieltjes transform) and Theorem 4.1 
(deterministic approximation of the mutual information of a wireless chan- 
nel). Theorems 2.4 and 2.5 are stated in Section 2. Section 3 is devoted 
to some examples and remarks while Section 4 is devoted to applications 
to digital communication; in particular. Theorem 4.1 is stated in Section 4. 
Proof of Theorem 2.4 is established in Section 5 while Theorem 2.5 is proved 
to the main extent in Section 6. Technical results and some of the proofs are 
collected in the Appendix. 
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2. The deterministic equivalent: main results. 

2.1. Notation, assumptions and preliminary results. Let N = N{n) be a 
sequence of integers such that hm„^oo = c E (0,oo). We denote by i the 
complex number and by Im(z) the imaginary part of z G C. Consider 
an X n random matrix Yn where the entries are given by 

where X^j and {aij{n)) are defined below. 

Assumption A-1. The random variables {XJ^; I <i < N ,1 < j <n,n> 
1) are real, independent and identically distributed. They are centered with 
E(Xj")^ = 1. Moreover, there exists e > such that 

E|X,'5f +^ < oo, 

where E denotes expectation. 

Remark 2.1. We can (and will) assume in the proofs that e < 4 without 
lack of generality. 

Assumption A-2. There exists a nonnegative finite real number <Tmax 
such that the family of real numbers ((Tjj(n), 1 < i < N ,1 < j < n,n > 1) 
satisfies 

sup max |aij(n)| < fJmax- 

n>l 1<«<JV 
l<j<n 

Denote by var(Z) the variance of the random variable Z. Since var(YJ") = 
afj{n)/n, the family {aij{n)) will be called a variance profile. 

Let An = {Afj) he an N X n real deterministic matrix. We introduce the 
N X n matrix 

For every matrix M, we will denote by M'^ (resp. M*) its transpose (resp. its 
Hermitian adjoint), by Tr(M) [resp. det(M)] its trace (resp. its determinant 
if M is square) and by F^'^^^ , the empirical distribution function of the 
eigenvalues of MM'^ . The matrix /„ will always refer to the n x n identity. 

Let A be an n X n matrix with complex entries. If A is Hermitian and 
positive semi-definite (i.e., z*Az > for every z G C"), we write A>0. li A 
and B are Hermitian matrices, A> B means A — B > 0. 
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We define the matrix Im(A) by 

Im(A) = i(A-^*). 

Let B be the Borel cr-field on M. An nx n matrix-valued measure M on B is 
a matrix- valued function on B, the entries of which are complex measures. 
A positive matrix-valued measure M is such that for every Ag B, M(A) > 
[i.e., M(A) is a Hermitian and positive semi-definite matrix]. In this case, 
the diagonal entries are nonnegative measures and for every z G C", z*Mz is 
a scalar nonnegative measure. For more details on matrix-valued measures, 
the reader is referred to the standard textbook [20] . 

In this article, we shall deal with several norms. Denote by || • || the Eu- 
clidean norm for vectors. In the case of matrices, the norm || • ||sp will refer 
to the spectral norm and if Z is a complex-valued random variable, denote 
by \\Z\\p = (E|Z|P)^/*' for p > 0. If A:" is a topological space, we endow it 
with its Borel cr-algebra B{X) and we denote by V{X) the set of probability 
measures over X . 

We denote by C" = {z G C,Im(z) < 0} and by C+ = {z G C,Im(z) > 0}. 

Assumption A-3. Denote by the fcth column of A„, by a" its ^th 
row and by 

amax,n = max( || a^ || , ||a^ || ; 1 < A; < n, 1 < £ < iV). 
We assume that 

^max = SUpaniax,n < OO. 
n>l 

We will frequently drop n and simply write aij, A, Y , a^, and so on. 
If (ai , . . . , afc) is a finite sequence of real numbers, we denote by diag(ai , . . . , 
Uk) the k xk diagonal matrix whose diagonal elements are the Oj's. Let 

{2.1) Dj = diag(a,^- {n)-l<i<N) and A = diag(CT^j (n) ; 1 < j < n) . 

1 /2 ~ 1 /2 

We will denote by D- = diag((Tij,i < N) and = diag(o"ij,j < n). Let 
be a probability measure over M. Its Stieltjes transform / is defined by 

m=f^, zee-. 

Jr X - z 

We list below the main properties of the Stieltjes transforms that will be 
needed in the sequel. 

Proposition 2.2. The following properties hold true: 
1. Let f be the Stieltjes transform of fi, then 
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- the function f is analytic over C~^, 

- ifzeC+ then f{z)eC+, 

- the function f satisfies: \ f{z)\ < j^^^ and Im(j^) < — Im(z) over , 

- if 00,0) = then its Stieltjes transform f is analytic over C — M"*". 

Moreover, z E C"*" implies zf{z) S C"*", and the following controls hold: 

r|Im(z)|-\ i/zGC-M, 
(2.2) \f{z)\<l\z\-\ ,/^e(_oo,0), 

[ (dist(z,]R+))-i, i/zGC-M+, 

where dist stands for the Euclidean distance. 

2. Conversely, let f he a function analytic over C"*" such that f{z) S C"*" 
if z ^ C"*". // liniy^+oo —iyf{iy) = 1, then f is the Stieltjes transform of a 
probability measure ^ and the following inversion formula holds: 

b 



/i([a,6])= lim - / Im/(e + ir/)dC, 



whenever a and h are continuity points of fi. If moreover zf{z) G for z 
in C"^ then, /i(M~) = 0. In particular, f is given by 



fiz) 



njdX) 
X-z' 



and has an analytic continuation on C — M"*". 

3. Let F be an n X n matrix-valued function over such that 

- the function F is analytic, 

- the matrices Im F{z) and Im zF{z) satisfy 

ImF(z)>0 and lmzF{z)>0 forz£C+. 

Then there exists an n x n matrix C > and a positive nx n matrix-valued 
measure /i over M"*" such that 

Fiz)=C+( ^ wUhTrf ^<oo. 
Jk+ X-z Jr+ 1 + a 

4. Let Fn and F be probability measures over M and denote by /„ and f 
their Stieltjes transforms. Then 

(yz G C+, fn{z) fiz)) ^ P„ ^ P, 

\ n—>oo 1 n—>oo 

where T> stands for convergence in distribution. 



As item 3 is perhaps less standard than the other items, we give a proof 
in Appendix A. These results can be found in several papers without proof 
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(see, e.g., [3], pages 64-65). 



In the sequel we shall denote by 5(M+) the set of Stieltjes transforms of 
probability measures over M^. Proposition 2.2 part 2 gives necessary and 
sufficient conditions for / to be in 

We list below some useful properties of matrices that will be of constant 
use later: 



Proposition 2.3. 1. Let B he an N ^ n matrix and let C he an n^k 
matrix, then 

||i?C||,p<p||sp||C||sp. 

2. Let B he an N xn matrix, let y and x he respectively iV x 1 and n x 1 
vectors, then 

ll^llsp = sup{|y^-Bx|, ||y|| = ||x|| = 1}. 



In particular. 



I L^ij I ^ 1 1 1 1 sp ■ 



3. If B and C are square and invertible matrices, then 

5-1 - = -B-^{B - C)C-^. 

If B and C are resolvents (see below), then this identity is known as the 
resolvent identity. 

4- Let X he an N X n matrix and denote hy Q{z) = {Qij{z)) the resolvent 
ofXX^, that is, Q{z) = {XX^ - zl)~^ . 

(a) We have ||(5(-2)||sp < |Im(2;)|-i and in particular \Qij{z)\ < \ \m.{z)\~^ . 

(b) If D is an N X N diagonal matrix with nonnegative diagonal elements, 
then 

Although very simple, the identity in Proposition 2.3 part 3 is extremely 
useful and will be of constant use in the sequel. 

2.2. The deterministic equivalent: existence and asymptotic behavior. Let 
us introduce the following resolvents: 

Qn{z) = (SnS^ - zIn)~^ = {qij{z))l<i,j<N , Z£C- M+, 

Qn{z) = (S^S„ - Zln)'^ = {qij{z))i<ij<n, Z G C - M+. 

The function fn{z) = j^Ti Qn{z) is the Stieltjes transform of the empirical 
distribution of the eigenvalues of The aim of this paper is to get 



(2.5) 
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some insight on /„ by the means of a deterministic equivalent. We wih 
often drop subscripts n. In Theorem 2.4, we prove the existence of matrix- 
valued deterministic functions T{z) and T{z) which satisfy a system of A^ + n 
functional equations. In Theorem 2.5, we prove that asymptotically, 

lTrQ(z)~i-Trr(z) and -Tr Q{z) ^ -Tifiz) 
Jy N n n 

in a sense to be defined. 

Theorem 2.4. Let An he an N x n deterministic matrix. The deter- 
ministic system of N + n equations, 

(2.3) Mz) = ~^ ~ ~ forl<i<N, 

z{l + {l/n)TTD,T{z)) 

(2.4) i^jiz) = for l<j<n, 

z{l + {l/n)TT DjT{z)) 

where 

^{z) = diag{'ilji{z),l<i<N), 
^{z) = diag{Tpj{z),l<j<n), 

T{z) = {^~\z) - zA^{z)A^)-\ 
f{z) = {^~\z) - zA'^^{z)Ay\ 

admits a unique solution (V'l, . . . , i^Nii'i, ■ ■ ■ , ^n) in 5(IR+)^"'"". 

Moreover, there exist a positive N x N matrix-valued measure fi = {mj) 
and a positive nxn matrix-valued measure fi = {jiij) such that 

I^l{M+)=In, /i(M+)=/„ and 

Tiz) = [ nz)=f ^ 

Jm.+ X - z Jm.+ a - z 

for z €C — M+. In particular, ■^TrT(z) and ^TrT{z) are Stieltjes trans- 
forms of probability measures. 

Proof of Theorem 2.4 is postponed to Section 5. It relies on an iteration 
scheme and on properties of matrix-valued Stieltjes transforms. 

Theorem 2.5. Assume that Assumptions A-1, A-2 and A-3 hold, then 
the following limits hold true almost everywhere: 

lim ( — TrQ(z) - — Trr(z) I =0 VzGC-M+, 
lim (-T):Q{z)--TYf{z)) =0 VzgC-]R+, 

n~^oo,N/n~*c\n U J 
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Remark 2.6 (Limiting behavior). There are two well-known cases where 
the empirical distribution of the eigenvalues of converges toward a 

limit expressed in terms of its Stieltjes transform: The case where the vari- 
ance profile (Tij = (T is a constant [5] and the case where the centering matrix 
A (which can be rectangular) has elements equal to zero outside the diag- 
onal [13] (with F"^""^" converging to a probability measure in both cases). 
Interestingly, one can obtain discretized versions of the limiting equations 
in [5] and [13] by combining (2.3)-(2.5). 

Corollary 2.7. Assume that Assumptions A-1, A-2 and A-3 hold. 
Denote by and 7r„ the probability measures whose Stieltjes transform are 
respectively j^TiQn and -^TrT^. Then the following limit holds true almost 
everywhere: 

POO roc 

/ /(A)P„(dA) - / /(A)7r„(dA) 0, 
Jo Jo 

where f : M"^ ^M. is a continuous and bounded function. The same results 
hold for the probability measures related to ^TicQn and ^TrT^. 

Proof of Corollary 2.7. We rely on the fact that (P„) is almost 
surely tight and that (vr^) is tight. In order to prove this, it is sufficient to 
prove 

sup^y AP„((iA)^ < oo a.s. and sup(^J A7r„((iA)^ < oo. 
One has 

J XFnidX) = ^TVS^S^ 

< ^TvAnAl + ^TrYnY^ 

(a) 9n-2 

— max 

where (a) follows from Assumption A-3. The almost sure tightness of (P„,) 
follows then from the fact that ^ij converges almost surely to one. 

The fact that sup„(/ A7r„((iA)) < oo is established in Lemma C.l. 

Let us now consider a realization for which the conclusion of Theorem 
2.5 holds, and (P^) is tight. Then for any subsequence n' of the natural 
numbers, there exists a further subsequence n" for which both P„// and 7r„" 
converge in distribution to, say P* and vr*. Then their Stieltjes transforms 
also converge and, from Theorem 2.5, to each other. Therefore P* = vr* and 
the conclusion holds on (n"). Since this is true for any subsequence (n'), the 
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conclusion of Corollary 2.7 holds for all n, for this realization, which occurs 
with probability one. □ 

We are indebted to the referee for the proof of Corollary 2.7 which is much 
simpler than the original one. 

Remark 2.8 (Concentration and martingale arguments). If one is inter- 
ested in proving that the empirical measure is close to its expectation, one 
can rely on concentration arguments (see, e.g., [11]) at least when the entries 
are Gaussian, bounded or satisfy the Poincare inequality. One can also rely 
on martingale arguments, regardless of the nature of the entries (as long as 
they are independent and satisfy some mild moment assumptions, see [10], 
Chapter 16, and also [6]). The purpose here is to provide a "computable" 
deterministic equivalent which is not expressed in terms of expectations. 
In fact, although expectations can be computed by Monte Carlo methods, 
these methods quickly imply a huge amount of computations when the size 
of the matrix models increases. 

2.3. Outline of the proof of Theorem 2.5. The proof relies on the intro- 
duction of intermediate quantities which are the stochastic counterparts of 
^, T and r. Denote by 



(2.6) 

(2.7) 
and by 
(2.8) 



hiiz) = forl<i<iV, 

z{l + {l/n)TrDiQ{z)) " " 

B{z) = d\ag{hi{z),l<i<N), 

'^^'^ = z{l + {l/n)\rD,Q{z)) '''' ^ ^ ^' ^ 
S(z)=diag(6,(z),l<j<n), 



R{z) = iB-\z) - zAB{z)A'^y^ and 
R{z) = {B-\z) - zA^B{z)Ay\ 



The introduction of the quantities R, R, T and T can be traced back to the 
work of Girko. We first prove that j^Tv Q{z) ~ j^TrR^z) and jjTtQ{z) ~ 
■^Tr R{z) in Lemma 6.1. These computations, quite involved, are along the 
same lines as the computations by Dozier and Silverstein in [5]. We then 
prove that jfTrR(z) ~ TrT(2:) and Tr R{z) ~ Tr T'(2:) in Lemma 6.6. 
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3. Examples and remarks. 

3.1. An example where the spectral measure of T,nT,'^ does not converge. 
Let Yji he a 2n X 2n matrix such that 

'Wn 0' 





Yn 



where Wn is n x ?i and Wljj = the Xjj's being i.i.d. Matrix Yn can be 
interpreted as a matrix with a variance profile a'^{x,y) in the sense that 

^'■'^ where a(.,,) = (i' if ^ [0, 1/2] x [0, 1/2], 

^ \0, else. 

Consider now the following 2n x 2n deterministic matrices Bn and Bn de- 
fined as 

Then both F^"^^ and F^"^" converge toward ^6o{dx) + ^Si{dx). It is well 

known that the limiting distribution of F^"^" and F^^"'^^"^^^"'^^"^'^ exist 
and are nonrandom probability distributions. The first limiting probability 
distribution is known as Marcenko-Pastur's distribution, and we denote it 
by IPj;;ip. Denote by Pcub the limiting distribution (This probability distribu- 
tion is sometimes referred to as the "cubic law" since the limiting Stieltjes 
transform of the normalized trace of the eigenvalues of {Wn + In){Wn + In)'^ 
satisfies a cubic equation (see [10], Chapter 2 for details), hence the notation 
Pcub.) of (for details, see [9, 13, 19, 25]). The following 

holds true: 



Proposition 3.1. Let T„ = (y„ + Bn) and let T„ = {Yn + Bn). Then 

I Pcub + 



n— >oo 



n— >oo ^ ^ 



Consider in particular the 2n x 2n random matrix defined by 

^('^n, ifn = 2p, 
\fn, ifn = 2p+l, 

that is, Tin = (Yn + An) where An is alternatively equal to Bn and to Bn, 
then F'^"'^" does not admit a limiting distribution. 
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Proof. Let {Wn + /„)(W„ + Inf = Un^nU^ where A„ = diag(Ai, 1 < 
i <n) and Un is orthogonal. Then ^ ^A; ^cuh • Now 

TT-(Un \ /A„ 0\ (Ul \ 
'"'"-1^0 /J 1^0 ojlvO /J' 

and F^"^" = ^(X]r=i '^A^ + '^'^o) which converges toward ^Pcub + ^<^o- One 
can prove similarly the second statement of the proposition. The absence of 
convergence for F^"^^ follows. □ 

Remark 3.2. In view of Proposition 3.1, F^"^" does not have a limit 
since it oscillates between two different limiting distributions, despite the 
fact that: 

(i) the variance profile is the sampling of a given function, 

(ii) F"^"^^ converges toward ^6o + ^Si. 

The general model S„ = (y„ + An) where is diagonal has been studied in 
[13]. In particular, an assumption is required in order to insure the conver- 
gence of F^"^" . In the setting of the present example [where Yn is defined 
by (3.1) and A„ is alternatively equal to Bn or this assumption writes: 
Denote by {an, 1 <i < 2n) the diagonal elements of An, then if the empirical 
(deterministic) distribution 

2n 

2 = 1 

converges toward a probability distribution with a compact support, then 
F " converges toward a deterministic distribution which is properly de- 
scribed via its Stieltjes transform. 

This assumption expresses a "coupled convergence" between the variance 
profile and the diagonal elements of An = diag(oji) and one can easily check 
that it is not fulfilled here. In fact, denote by {bn, 1 <i < 2n) the diagonal 
elements of matrix Bn and by {bn, 1 <i < 2n) those of matrix Bn- Then if 
n = 2p, 

^ 2n 

^^ = ^J2 ^{i/2n,bl) 1 [0,1/2] {x) dx (g) Si{dy) + 1 [1/2,1] {x) dx ® 6o{dy), 

i=l 

while if n = 2p + 1 , we have 
^ 2n 

= \i/2n,bl) ^"i;;^ l[o,i/2] {x) dx ® do{dy) + l[i/2,i] {x) dx Si{dy). 

4 = 1 

This, in particular, does not allow the convergence of Hn- We believe that 
this is the main reason for which does not converge in the present 

example. 
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3.2. The case of a separable variance profile. In this section, we consider 
the case where the variance profile is separable, that is, 

cr| = didj, l<i<N,l<j <n. 

Note that this occurs quite frequently in the context of MIMO wireless 
systems. The system of + n equations defining ^{z) and ^{z) takes a 
much simpler form: In fact it reduces to a pair of equations. To see this, 
we denote by D and D the diagonal matrices D = diag{di, ... ,dj\f) and 
D = diag((ii, . . . ,c?„). Note that Dj = djD and Di = diD. Denote by 5{z) 
and 6{z) the functions defined on C — M"*" by 



(3.2) 



5{z) = -lk{DT{z)), 
n 

~5iz) = -TviDTiz)). 
n 



It is obvious that ^{z) = + 5D)~i and ^{z) = + 5D)-^. There- 
fore, the solutions of the system of + n equations (2.5) can be expressed 
in terms of functions 6{z) and 6{z). Using the fact that T{z) = (^'(z)~^ — 
zA^{z)A'^)-^ and f{z) = {^{z)"^ - zA^^{z)A)-^ and the definition (3.2), 
we get that 5{z) and 5{z) are solutions of the following system of two equa- 
tions: 

6{z) = - Tt[D{-z{I + D5) + A{I + D6)~^A'^y^], 
n 

6(z) = - Tt[D(-z(I + D5) + A^(I + D6r^Ar^]. 
n 

Of course. Theorem 2.4 implies that this system has a unique solution in 
the class of functions {6,6) where 6 and 6 are Stieltjes transforms of positive 
measures (/i,/2) on M+ such that ^(M+) = ^Tr(D) and /i(M+) = ^Tr(D). 
Moreover at a given point z = —cr'^, it can be shown that the following 
system: 



(3.3) 



(3.4) 



- TT\D(a^(I + Dk) + A(I + Dk)-^A^)-\ 
n 

-Ti\D{(7'^(I + DK)+A^{I + DRY^Ay\ 
n 



has a unique pair of solutions (k, k), k > 0, k > 0. This, of course, is of 
practical interest in order to compute the approximant (4.3). 

Remark 3.3 (Limiting behavior, followed). The case with a separable 
variance profile illustrates quite well the kind of assumptions one needs in 
order to obtain a limiting behavior for the probability distribution of the 
eigenvalues of Sn^^n • If Tr (5(2:) converges then so does TrT(z) by The- 
orem 2.5. This in turn should imply the convergence of the right-hand side 
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of (3.3). But, unless D = I and D = I ox A pseudo-diagonal, such a con- 
vergence relies on an intricate assumption between the limiting distribution 
of D, D and the spectral properties of A. Otherwise stated, it is far from 
obvious to obtain a limiting behavior in the noncentered case even when the 
variance profile is separable. 

3.3. Additional information about the bias of the deterministic equivalent 
in a simple case. Let X„ be an x n matrix with i.i.d. entries and de- 
note by A„ = diag(Ai, . . . , Xjy) where jf J2f=i ^x'^ ~^ H{dX). We consider the 
matrix 

In 



Assume that ^ ^ c G (0, 1) then it is well known (see, e.g., [25]) that the 
spectral distribution of YnXn converges toward a probability distribution /x 
whose Stieltjes transform / satisfies 

We shall prove, in this simple setting, that the deterministic approximation 
/n(-z) = Z^i^i ''/'i(^) of the empirical Stieltjes transform mn{z) = Tr(y„y^ — 
zl)~^ satisfies the following discretized version of (3.5), where H{dX) is re- 
placed by Hn{dX) = i Y.f=i {dX) and c, by c„ = ^: 

1 ^ 1 

In this simple case, the variance profile cr(i,j) is equal to Aj and from equa- 
tions (2.3) and (2.4), we can write: 

''''^'^^ z{l + {Xj/n)EU^j{z)y "^'^'^^ z{l + {l/n)EliXjUz)) 
or equivalently 

V.(^)(i + f Ev;.(^)) =-K i^j{z)(^i + l^j:x!Mz)^=-l. 

Eliminating (after immediate manipulation) J2i ^ii^i J2j '4'j both equa- 
tions yields 

N 



1 / 1 \ 

where c„ = ^. The deterministic approximation fn{z) thus satisfies (3.6). 



16 W. HACHEM, P. LOUBATON AND J. NAJIM 

One can notice that the centering element in the CLT estabhshed by Bai 
and Silverstein in [2] is based on this kind of deterministic approximation. 
Moreover, the hmiting bias 

l3n = Emn{z) - fn{z) 

is computed in ([2], Section 4) when n — > oo. According to these results, it 
turns out that: 

1. The bias /?„ is of order if the entries of Xn are real, ^Xf^ = 3. 

2. The bias [3n is of order o(^) if the entries of Xn are complex, E|Xjj|^ = 2 
and EX2. = 0. 

3.4. The complex case and the case of a Gaussian stationary field. 

The complex case. In the case where the entries of X^ and are com- 
plex, the interest lies in the spectrum of (Yn + -^njiXn + ^n)* • All the results 
established in the real case hold in the complex case, with the following slight 
adaptations: 

1. In Assumption A-1, the random variable Xf^ is complex. Replace "^Xf^ = 
1 by ElATijf = 1 and EX^- = 0. 

2. In Assumption A-3, the entries of are complex. 

3. In Theorem 2.4, T and T must be replaced by the following: 

T{z) = {^~^{z) - f{z) = {^^^{z) - zA*^{z)A)-^. 

(3.7) 

4. In Theorem 2.5, Qn and Qn denote 

Qn{z) = (SnS* - zInY^, Qn{z) = (S;^Sn " zln)~^ ■ 

Similar adaptations must be made in the remainder of the article. 

The case of a Gaussian stationary field. Following [12], it is possible to 
rely on Theorem 2.5 to provide a deterministic equivalent in the case where 
matrix is replaced by Zn = (-^^u^) where 

(3.8) ZIj2 = ^ E Hh,k2)U{ji-h,j2-k2), 

where h is a deterministic complex summable sequence and {U{ji,j2);{ji, 
J2) G Z^) is a sequence of independent complex Gaussian random variables 
with zero mean and unit variance, that is, EC/ (71,^2)^ = 0, E|C/(ji,j2)P = 1 
[otherwise stated, C/(ii,j2) = V^ + iVF, where V and W are independent and 
AA(0, ;^)-distributed]. Denote by 

$(ti,t2)= ^ /i(4,£2)e2-^('i*i-'2*2). 
(£1/2)6^2 
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Function |<I>| is in particular real and bounded due to the fact that /i is a 
summable sequence. Denote by Fp the pxp Fourier matrix Fp = {F^^ j2)o<jij2<p 
defined by 



We can now state: 

Theorem 3.4. Let S.„ = Z„ + Bn where Z^s entries are given by (3.8) 
and Bn is a complex N x n matrix which satisfies Assumption A-3. Denote 
by Q{z) = (SnS* - zInT^ and Qn{z) = (S-^^Sn - zln)~^ ■ Then 

hm (lTrQ(z)-lTrr(z))=0 

and 

lim (—TvQ{z) TtT{z] 

n^oc,N/n-^c\n n 

for every z G C — where T and T are given by (3.7), where the variance 
profile aij{n) must be replaced by \^{jf,n)\, and matrix A by F^BF*. The 
convergence holds in probability. 

Lemma 3.3 in [12] [whose assumption is fulfilled due to Assumption A-3] 
together with Proposition 2.2 part 4 immediately yield the desired result. 

4. Application to digital communication: approximation of the mutual 
information of a wireless MIMO channel. The mutual information is the 
maximum number of bits per second per Hertz per antenna that can be 
transmitted reliably on a MIMO channel with channel matrix H^. It is 
equal to 

(4.1) Cn{a^) = -^Elogdet(^/^ + 

where cr^ represents the variance of an additive noise corrupting the received 
signals. The mutual information C„((J^) is related to ■^Tr(i/„ff* + (T^/tv)""^ 
by the formula 

^ = lETV(i/.if: + .^/^)--J, 

or equivalently by 

Cn{T^) = j ^ (^- - -ETr(i?„i7: +c^/^)-ij du:, 
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which fohows from (4.1) by Fubini's theorem. In fact, ^ — j^Tr{HnH* + 

ujI)^^ = — Jq° '^A+l^^ — ^' '^here stands for the empirical distribution 
of the eigenvalues of HnH*. In certain cases, channel matrix Hn is unitar- 
ily equivalent to S„ = 1^ + An which has complex entries. Without loss 
of generality, we shall work with matrix S„ with real entries in order to 
remain consistent with the general exposition. Showing that ■^Tr(SnS^ + 
ujIn)~^ — ;^Trr„(— w) for the deterministic matrix-valued function Tn{z) 
defined in Theorem 2.4 allows one to approximate Cn{o''^) by Cn{o''^) = 
/(t2°°(^ ~ j^'^Tn{—iji>))duj. This approximant can be written more explic- 
itly. 



Theorem 4.1. Assume that Assumptions A-1, A-2 and A-3 hold anc 
denote by 

r+oo /II \ 

(4.2) Cn{a^)= ---Trr„(-a;)Ua;, 



.OJ N 

where T is given by Theorem 2.4. Then the following limit holds true: 

n^+oo,N/n^c 

where cr^ E M"*" . Moreover, 



C„(a2) = llogdet 



2 +A^{-a^)A^ 



(4.3) 

+ 1 logdet _ ^ 

Proof of Theorem 4.1 is postponed to Appendix C. In certain cases, the 
study of the behavior of Cn{o''^) is simpler than the behavior of Cn{o''^), and 
allows one to get some insight on the behavior of the mutual information of 
certain MIMO wireless channels (see, e.g., [7] for preliminary results). 

Remark 4.2. Equation (4.3) has already been established in the zero 
mean case {An = 0): the centered case with no variance profile has been 
studied by Verdii and Shamai [27] , the centered case with a variance profile 
by Sengupta and Mitra [22] (with a separable variance profile) and Tulino 
and Verdu ([26], Theorem 2.44). 



5. Proof of Theorem 2.4. We will first prove the uniqueness of the solu- 
tions to the system (2.3) and (2.4) within the class 5(M"'") and then prove 
the existence of such solutions. The following proposition will be useful. 
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Proposition 5.1 (Matrix- valued Stieltjes transforms). Let An be an 
N xn real matrix. Let ji E 5(M+), 1 <i< N, and 7j G 5(M+), l<j<n. 
Denote by 

r = diag(7i, 1 < j < A^), f = diag(7,-, l<i<n), 
T = (r-i - zATA^)-^, T = (f-^ - zA^TA)-^. 

Then: 

1. The matrix-valued functions T and T are analytic over C — M"^. 

2. If z £ C~^, then ImT(2;) > and lmzT{z) > 0. Moreover, there exists 
a positive N x N matrix-valued measure fi = {fJ,ij) and a positive n x n 
matrix-valued measure p, = {fiij) such that 

/[/(R+) = In, /i(M+) = /„ and 
T(z) = / T(.) - t 



A — z ' 7iR+ \ — z 

for ^; E C — M"^. In particular, jj Tr T{z) and ^ Tr T(z) belong to S{ 
3. The following inequalities hold true: 

VzGC+ T{z)T*{z) < and t(z)t*(z) < 



(Imz)2 (Imz)^ 
4. Lei Dj and Dj be defined by (2.1). Denote by 



z(l + (l/n)Tr(AT(z)))' 



z(l + (l/n)T\-(D,T(z)))^ 



where 1 < i < N and 1 < j < n. Then and 7^ are analytic over 
C - M+ and belong to 5(M+). 

Proof. We only prove the stated results for T, the adaptation to T 
being straightforward. 

Since ^i{z) S 5(R^), we have '^i{z) ^ over C — M"*". In particular, F"^ — 
zATA^ is analytic over C — M^. In order to prove that T is analytic, it 
is sufficient to prove that det(r-i - zATA^) / for every z E C — M"*" or 
equivalently, to prove that if /i G is such that (F^^ - zArA'^)h = 0, then 
h must be equal to zero. Let h be such that (F~^ — zATA'^)h = 0, then 

( h*{T-^ - zAtA^)h = 0, 
^ • ' \h*{T-^* - z*Ar*A'^)h = 0. 
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We study separately the case z G (—00, 0) and the cases z G C"*" and z G C~ . 
If z G (— c«,0), then r^^(z) > \z\lfq and —zt(^z^ > and the first equation 
in (5.1) yields: 

= h*{r~\z) - zAr{z)A^)h > \\hf\z\ 

which implies /i = 0. If z G C"*", then Im(zr(z)) > and Imr^^(z) < — Im(z)I„ 
by Proposition 2.2 part 1. The system (5.1) yields h* (lm.{T~^ (z)) — AIm.{zr{z)) x 
A'^)h = 0. Necessarily, 

= h*{lm{T^\z)) - Alm{zt{z))A'^)h < -{lmz)\\hf, 

which yields h = 0. If z £ C~, the same kind of arguments yields h = 0. 
Therefore, T is analytic over C — M"*" and part 1 is proved. 

In order to apply Proposition 2.2 part 3, we now check that ImT(z) > 
and ImzT(z) > for z G C"*": 

^ ^ 2i 

= T{z){-lmT-\z) + Alm{zr{z))A'^)T*{z) > 0, 

where (a) follows from the resolvent identity and the last inequality follows 
from the fact that 7^ and 7j belong to 5(M~'~) and from Proposition 2.2 
part 1. 

One can check similarly that Im(2;T(z)) > 0. Therefore Proposition 2.2 
part 3 yields T{z) = C + ^x^' '^here C = iCij) \s an N >i N nonnegative 
matrix and ^ \s an N x N matrix-valued measure. Let us now prove that 

(5.2) lirn iyT(iy) = -/^. 

J/— >+oo 

First write 

(5.3) T(2) = (r-i(z) - zAt{z)A^y^ = {In - T{z)Azt{z)A^)-^T{z). 
Since F and F are diagonal matrices whose entries belong to 5(]R"*"), we have 

(5.4) lirn iyF(iy) = -/at, liin iyf(iy) = 

y— >+oo J/— >+oo 

Therefore, it is straightforward to show 

(5.5) ||F(i2/)^(iyf(iy))^^||sp — > 0. 

Plugging (5.4) and (5.5) into (5.3), we get (5.2). 

In order to prove that C = and ^(M+) = In, we introduce the complex 
functions: fi{z) = Tjj(z),l <i < N. These functions are analytic over C"*" 
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and satisfy: Vz G C"*", Im/j(z) > 0,lmz/j(z) > 0. Moreover (5.2) implies that 
limj^_>+oo —iyfi{iy) = 1; therefore Proposition 2.2 part 2 yields the existence 
and uniqueness of Ui G 'P(IR"*") such that 

On the other hand, 



A 



= Tu{z) = Cu+ f f^f^, l<z<N. 
Jr+ \ — z 

From this, we deduce that /ijj = z/j G ^'(M"^) and that Cu = 0. This im- 
plies that C = since C > 0. As ^ is a positive matrix-valued measure, 
l/ifcfl < ^(m +/^«), =0 and |/ifc£|(M+) < cx). Assume that ki^^. 

Equation (5.2) yields 

(5.6) lirn iyTfc£(iy) = 0. 



But 



f f 

iyTke{iy) = - / t^- — ^m{dX) + i / -r^ — KfikddX). 



A^ + y2 ^2 _j_ y2 

2 

In particular, (5.6) implies that liuiy^^oo 'i l^keidX) = 0. But the 

2 

dominated convergence theorem implies that limy^+oo/K+ x^\y^ fJ-keidX) = 
/ifc^(M+). Necessarily fikeO^'^) = 0, /^(M"*") = In and part 2 is proved for 
z £ C^. Since both T and / are analytic over C — M"*^ and coincide 
over C"*", they are equal. Taking the trace, we immediately obtain that 
;^TrT(z) G 5(M"'') and part 2 is proved. 

In order to prove part 3, we shall prove the two following inequalities: 

2i Imlz) 
T(.) -^^(^r >T(,)T*(z)(Im(.)). 

The first one is straightforward: 

T(z)-T(z)* f fiidX) 



ImT(z) = Im(2;) / 



2i Jr+ |A - ' 

Since fJ-{A) > for every Borel set A, so is / /(A)/x((iA) whenever / > 0. 

1 _ 1 



Applying this property to /(A) = . — jv^-ri > and using the fact that 



//(M"*") = In, we obtain 



|A — ~ Im(z) 
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which yields < jj^^. Using the resolvent identity, we get 



2i ^^^i '-^ + Alm{zT{z))A' ]T*{z). 

Let z G C"*", then 



2i °V \li{z) 

by Proposition 2.2 part 1. This and the fact that lm{zT{z)) > yields the 
second inequality in (5.7) and the proof of part 3 is complete. 

Concerning part 4, we only prove the statement for G 5(M''') (since 

the same arguments yield the desired results for 7] )• For this we rely on 
the characterization given in Proposition 2.2 part 2. The quantity z{l + 
^ Tr(l)jT(2))) is analytic over C"*", and its imaginary part is greater than 

Im(z) there. Therefore, it does not vanish over C~^, and 7^^ is analytic over 

C+. Similar arguments yield that 7^^ is analytic over C~ and (— C)O,0). Let 

us now compute Im7) {z): 

^f\z)-^f\zr Im(z) 



2i |z(l + (l/n)TV(AT(z)))|2 

Th:[b,{zT{z) - z*T{zr)/{2i)] ^ ^ 
n|2(l + (l/n)Tr(AT(z)))|2 " 

for z G C+ by Proposition 5.1 part 2. Similarly, one can prove that Imz7j^^^ {z) > 
for z G C+. Finally, 

l + (l/n)Tr(AT(iy))2^-^°° 

since | Tr(AT(iy))| < nal^^J^T {\y)\\sp and ||T(i?/)||sp < i. The proof of part 4 
is completed and Proposition 5.1 is proved. □ 

Proposition 5.2 (Uniqueness of solutions). Assume that {ipi,ipj]\ < 
i<N,l<j<n)£ cS(R+)^+" and {(jiiAj] I <i < N ,1 < j <n) e 5(M+)^+" 
are two sets of solutions to the system (2.3) and (2.4). Then these solutions 
are equal. 

Proof. Denote by 

^{z) = d\ag{Uz), l<i<N), ^z) = diag((^,(z), 1 < i < n), 

T^{z) = {^-^{z) - zA^{z)A^)-^, f^{z) = (4>~H^) - zA^^z)A)-^, 
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^iz) = diagiiPiiz),l<i<N), 
T^{z) = i^^-^z) - zA^{z)A'^)-\ 
Let us compute tpi — (pi for z G C"*". 

ipi{z) - (pi{z) 



§(z) = diag('0i(z),l <i <n), 
r^'(^) = - zA^^{z)A)-\ 



+ 



z{l + (1/n) Tr z{l + {l/n)Tr DiT'l'{z)) 



If z G C+ then 



z(l + -TrDif^{z) 
n 



>lmUU+^i:TDif'l'{z)\ \ > Im( 



where (a) follows from Proposition 5.1 part 2. The same argument yields a 
similar bound for 4>i{z), that is |(^^^(z)| > Im(z), and 



(5.8) 



mz)-uz)\ < 



-T\-A(T'^-r^) 



n 



(Imz)^ 

A standard computation involving the resolvent identity yields 
TrA(r"^-T'^) 

= - Tr AT'^C'^"^ - ^"^ - zA^{^ - ^)A)f'f' 

= -TTDif^{^~^ - §"i)r'^ + zTr Ar'^(^^(^ - ^)A)T'^ 



N 



= Tr AT'^^^'(«' - ^)^~'f^' + ^ zTr AT'^a^afc(<^fc - Vfc)^ 

k=l 

(recall that is the kth row of ^4). We introduce the following maximum: 

M = sup{|0i - i;i\,\4>j -i>j\,l<i<N,l<j <n}. 
We also define dn by d„ = max(-^, ^). Note that (i„, > 1. Then, 
1 

n 

First notice that Proposition 5.1 part 3 yields ||T'^||sp < (same inequality 
for T^, and f^). Since ^ = -z{l + ^Tr Ar^(z)), we have 



■TrA(T*-r'^) <Md„aLxr*llsp||T^1|sp{||l>^1sp||^-isp + k|aLx,n}- 



||^"isp < |Z|(1 +dn<axr^Wllsp) < k| 1 + ^^n 



Imz 



2 , 
max 



One can show similarly that max( II <I> ^HspJI^' ^||sp,||^ ^||sp) < kKl + dn^j^ 

^yrt n v dn. I Z 



Therefore 
1 



(5.9) 



n 



< 



(Imz) 



z\ \ l + d„ 



Imz 
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Plugging (5.9) into (5.8), we obtain 

\Mz)-Mz)\<£{z)M{z), 



(5.10) 

where £{z) = '^''^"^^^ \\z\{l + ] + 



|2 J 2 f ^ 0-2x2 

(Imz)4 V^l^v "Im^; ' 
The same kind of computation yields 
(5.11) \^Pi{z)-Mz)\<£{z)M{z). 
Finally, gathering (5.10) and (5.11) we obtain 

M{z) <£{z)M{z). 

Let z be such that \ < 2. For Imz large enough, one has £{z) < ^, which 
implies that M(z) = 0. The analyticity of (j)i,Tpi,(f)j^tjjj implies that = tpi 
for 1 < i < and (f)j = ipj for 1 < j < n on C — M"*". Proposition 5.2 is proved. 
□ 

Proposition 5.3 (Existence of solutions). There exists 1 <i < 

N,l<j<n)£ cS(M+)^+" satisfying (2.3) and (2.4). 

Proof. We construct the desired solution by induction. Let 

il^P (z) = (-2) = for 1 < i < A^, 1 < i < n. 

Then ^pf^ and i^f^ belong to 5(M+). For p > 0, let 

(z) = ~^ , , for 1< i < iV, 

' z(l + (l/n)Tr(AT(p)(z))) 

-i^l'^"''"^^ (z) = — — — , , for I < j <n, 



where 



z{l + {l/n)TT{DjTiP){z))) 

^^P\z) = diag(V'l^^(z), l<i<N), 

¥p\z) = diag(Vif ^(z), 1 < i < n), 

T(p){z) = (^(P)-i(z) - zA^{z)^P^A^)-\ 

f(p){z) = {¥p^-\z) - zA^^^P\z)Ay\ 

By Proposition 5.1 part 4, (^^^^-i^j^^; l<i<N,l<j<n) are analytic over 
C - M+ and belong to cS(M+) for ah p>0. Denote by 

M^P^ = sup{ - V'?) I , I V^f -i^f\,l<i<N,l<j< n}. 
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The same computations as in the proof of Proposition 5.2 yield 

M(p+i)(z) <£:(z)M(p)(z), 

where S{z) = ^^^^{kl (l + d/^)\ aLx,n 

Let z G C"*" be such that | jjf^l < 2. For Imz large enough, ijj^f^ {z) and ipj^\z) 

are Cauchy sequences. Denote by tpi{z) and ipj{z) the corresponding limits. 

On the other hand, {'4^^^)p is a normal family over C — (see, e.g., [21]) 

since ip^^ is bounded on every compact set included in C — M"*" uniformly 
in p [see the third inequality of (2.2)]. Therefore one can extract, by the 
normal family theorem, a converging subsequence whose limit is analytic 
over C — M"*". Since the limit of any converging subsequence is equal to il^i{z) 

in the domain {z G C^, |2;|/| Imzl < 2, Imz large enough}, ij^'f^z) converges 
toward an analytic function on C — M"*" [that we still denote by 'il)i{z)\. One 

can apply the same arguments for ipj''\z). 

We now prove that V'j and ipj satisfy (2.3) and (2.4) where ^, ^, T and T 
are well defined. Let "^{z) = diag('0i(2;), 1 < i < N) and "^{z) = diag(^i(z), 
1 < i < ■"-) • The convergence of {ipl^^ )p immediately yields that 



lmtp'f\z) >O^ImV'j(^) >0, 

lmzip^\z) > ^ lmz^i{z) > 0, 
1 1 



|V?^(^)|<:^^ 1^.(^)1 < 



Im z Im z 

on C^. It remains to prove that lim^^oo —'wi^iiw) = 1- The same arguments 
as in the proof of Proposition 5.1 yield that both 

T{z) = {^'\z)- zA^{z)A^y'^ and f{z) = {^-\z)-zA^^{z)Ay^ 
are well defined on C — M"*". Moreover, 

T^P\z) T(z) and f^P\z) f(z) 

on C — M'^ . We can deduce that 

ipi(z) = forl<i<iV, 

z(l + (l/n)TV(Ar(z))) - - 



i/jjiz) = — — for 1 <j <n 



-1 

z{l + {l/n)T^{D,T{z))) 
on C"*", and hence on C — R^. Therefore, T^p\z)T^p\z)* < ^^^^yi implies 
that T{z)T{z)* < y2 and one can easily prove that 



lim -iyVi(iy) = 1 
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using the previous representation of ipi. Hence, ipi belongs to 5(M"'^) by 
Proposition 2.2 part 2. We can prove similarly that ipj £ 5(M^). Proposition 
5.3 is proved. □ 

The proof of Theorem 2.4 immediately follows from Propositions 5.1, 5.2 
and 5.3. 



6. Proof of Theorem 2.5. For the reader's ease, we recall some of the 
notation previously introduced. 
The resolvents: 

Q{z) = (SS^ - zl)-^, Q{z) = {Y^^T. - zl)~^. 

The stochastic intermediate quantities: 

B{z) = d\&g{bi{z),l<i<N), 



bi{z) 
b,{z) 



-1 



z{l + {l/n)TT{D,Q{z)))' 
-1 



z{l + il/n)TiiD,Q{z)))' 
R{z) = {B-\z) - zAB{z)A'^y\ 
The deterministic quantities: 

"^'^^^ " z(l + (l/n)Tr(AT(z)))' 
^{z) = diag{il;i{z),l<i<N), 
T{z) = - zA^{z)A'^)-\ 



B{z) = diag{bj{z),l<j<n), 
R{z) = {B~\z) - zA'^B{z)Ay\ 

^^■^^^ = z(l + (l/n)TrP,r(z)))' 
^{z)=diag{iPj{z),l<j<n), 
f{z) = - zA^^{z)A)~^. 



6.1. Evaluation of the differences j^Ti Q{z) — j^TrR{z) and ^Tr Q{z) — 
^TiRiz). 

Lemma 6.1. Let f/„ = diag(n", I <i < N) be a sequence of N x N di- 
agonal matrices and Un = diag(n", 1 < « < n), a sequence of n x n diagonal 
matrices. Assume that there exists Ki and Ki such that 

sup max < Ki < oo and sup max < Ki < oo. 



Then, if z £ C"*" , 



E 



(6.1) 



-TTiQ{z)-R{z))U 



l+e/4 



and 



E 



1 



n 



Tt{Q{z) - R{z))U 



n 



l+e/4 ■ 



n 
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In particular, for each z G C^, 

^T^iQiz)- Riz))U and -TiiQiz) - R{z))U ^ Q 

N n 

almost surely as N ^ oo and N/n ^ c > 0. 



The computations in the proof of Lemma 6.1, along the same hnes as 
the computations in [5], require some adaptation due to the fact that Yn 
has a variance profile. In this section, we state two intermediate results in 
Proposition 6.3 and Proposition 6.4. A sketch of the proof of Proposition 6.4 
is given while the full proof is postponed to Appendix B. 

We first need to introduce related matrix quantities when one column/row 
is removed. Denote by y^, aj and $,j the jth column of Y, A and S respec- 
tively, and by A^^^ and 'S^^^ the x (n — 1) matrices that remain after 
deleting the jth column from matrices A and S, respectively. Also, let 

Q(J)(2) = (S(J')s(-'')^-^/)"^ and g(^')(z) = (S(j')^S(j)-z/)-^ 

Let D^^^ be the (n - 1) X (n — 1) matrix where column j and row j have 
been removed: 



DY' = dia.g{al„l<£<n,£^j). 



Let 
(6.2) 

(6.3) 



R^^ 



z{l + {l/n)Ti{D\^^Q^i))) 

diag(6p'^(z),l <i< Af), 
-1 

z(l + (l/n)T¥(L>^QO-)(z))) 

diag(6?^(z),l<£<n,£/j), 
(S(j')"'(z) -zAW5(J')(z)^(J')^)~\ 
(S(j')"'(z) -z^(^')^5(^')(z)A(J'))-\ 



1 < i < A^, 



l<^<n,^/j, 



Note that any random matrix with superscript (j) is independent of y^. This 
fact will be frequently used in the sequel. By the matrix inversion lemma 
(see [16], Appendix A, see also [15], Section 0.7.4), we have 



(6.4) 



Q = qO) - -Lg(i)^^.^jQO) where aj = l + ^JQ'-^^^j ■ 
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Recall that = slj + yj . Following [5] , we introduce the following random 
variables: 



We note that aj = 1 + ujj + 2f3j + pj. 

Remark 6.2 (Some qualitative facts). We can roughly split the previous 
random variables into three classes: 

— The vanishing terms: The terms Pj, Pj and jj vanish whenever n goes to 
infinity; the main reason for this to hold true follows from Ey^ = (this 
is formally proved in Proposition B.l). 

- The quadratic forms based on y^: The behavior of the terms ujj and ujj 
can be deduced from the following well-known result: 

x'^Axr^-TiA, where x'^ = ^(Xi,...,Xn), 
n \ n 



the Xi being i.i.d., as long as x and A are independent [see (B.l) for the 
full statement]. Concerning a)j, R is not independent of yj however we 
shall see in the course of the proof that R can be replaced by which 
is independent of yj . Hence the previous result can also be applied to Coj . 
— The mysterious terms: Nothing is known about the asymptotic behavior 
of the terms pj and pj, except that these terms are bounded. Fortunately, 
these terms are always combined with vanishing terms, as we shall see in 
the sequel. 

Proposition 6.3. The following expression holds true: 
^ Tr(Q -R)U = Zl' + + Z^ + Z^ + Z^, 

where 

1 n ~ 
yn _ ^ Sr^ Ij 

^ B- 
j=l j 
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j=i J 
1 ^ 



Proof. Let us develop {R - Q)U: 
{R{z)-Q{z))U 

= Q{z){Q"\z) - R~\z))R{z)U 

= -zl- B-^{z) + zAB{z)A^)R{z)U 

= Q(^) I E ((1 + 4- + aj-yj + y^aj + y^yj) 



{B~\z)+zl) \r{z)U. 



In particular, 



1 11 
(6.5) - Tr(i?(z) - Q(z))C/ = - ^ Ty^- - _ Tr(i?-i + z/)i?(z)[/Q(z), 

i=i 

where 

Wj = {l + zbj)sijRUQaj + yj RUQaj +ajRUQyj + yjRUQyj. 

Using (6.4) and the fact that aj = 1 + ujj + 2/3j + Pj, a straightforward 
(though lengthy) computation yields: 



Wj = ^ + ^(1 - zhj{pj + ;3j)) + ^(1 + ^fei(l + /3j)) + 

Cfc 

Plugging this into (6.5), we obtain the desired result. □ 



Proposition 6.4. There exists a nonnegative constant K such that 

K 

K 



E|^3 < 



2+£/2 ■ 



n 



E|Zrr^/^<-A_, z = 2,4,5. 

We first give a sketch of the proof, based on the qualitative facts given in 
Remark 6.2. The full proof is given in Appendix B. 
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Remark 6.5 (Useful upper bounds). By Proposition 2.3 part 4(b), 6j, 
bg, h^p and V"p belong to 5(M"'") with probability one. Therefore, as a con- 
sequence of Proposition 5.1, 

with probability one. Furthermore, 

where (6.6) can be derived by applying the results in Sections 0.7.3 and 0.7.4 
in [15], for instance. Therefore, due to Proposition 2.3 part 4(a), we have 

z 

Imz 



(6.7) 



aAz) 



< 



Sketch of proof of Proposition 6.4. We loosely explain why the 
random variables (1 < i < 5) go to zero. We first look at Z". The term 7,- 
vanishes uniformly in the sense that one can prove that sup]^<j<„ WljWi ^ 

(cf. Proposition B.l). Since is bounded (see Remark 6.5), Minkowski's 

inequality implies that ||.Zf II4 = 0{^^). 

The terms Z2 and can be handled in the same way: The terms Pj 
and Pj vanish uniformly (cf. Proposition B.l). Furthermore, the terms pj, 
Pj, bj and aj^ are uniformly bounded. By consequence, Z2 and Z^ vanish 
for large n. 

The analysis of Z^ relies on arguments about quadratic forms. Denote by 
Xj = {Xi j, . . . , j)'^ . Then yj = -^Oy^-Kj and the quadratic form loj can 
be written: 

-. = yjQ(^)y. = ixJz.fQ(^)z)fx, 

As Xo and d]^'^Q^^^ d]^"^ are independent, and d]^'^Q^^^ d]^"^ is bounded, 
supi<j<„ \\ujj - i TrQDj\\2+s/2 = (see [1]). This implies that zbj{l + 

ujj) ~ —1 and that sup;^<j<,„ ||1 + 2:6^(1 -|-u;j)|| 2+^/2 = 0{-^). The term ctj^pj 
being bounded, Minkowski's inequality immediately yields 11^4 || 2+2/2 = ^^'^)' 
Let us now look at Z^. We have Uj = yjRUQ^^^yj. Perturbation arguments 
yield: 
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where "~" loosely means "asymptotically equivalent." Now approximation 
formulas for quadratic forms yield: 



Since {aj] 
(6.8) 



-zQjj, we have: 



1 \ - UJi 



1 



iV ^ a, 



n 



On the other hand, straightforward computation based on the mere defini- 
tion of bi yields: 



1 1/1 
(6.9) - TT{B-\z) + zlN)Riz)UQ{z) = --J2^zqjj{z)-TTDjRiz)UQ{z, 



N 



[see (B.17) for details]. Combining (6.8) and (6.9) allows to show that 

\\ZE\Ue/2 = Oi^). □ 

Proof of Lemma 6.1. Since ^Tr(Q(z) - R{z))U = + + + 
Z2 -\- Z^, Proposition 6.4 and Minkowski's inequality immediately yields 

2+e/2 



(6.10) 



E 



j^TviQ-R)U 



< 



n 



K2 
l+e/4 ■ 



Borel-Cantelli lemma implies that j^Tr{Q{z) — R{z))U almost surely. 
The assertion i Tt{Q{z) - R{z))U a.s. can be proved similarly. □ 

6.2. Evaluation of the differences j^TicR{z) — ■^TrT(2;) and ^Tr R{z) — 
^ TrT'(z). Recall that dn is defined by d„ = max(-^, ^).As^— >c, 0<c< 
+00, it is clear that sup^(i„ < d where d < +00. Denote by V the subset of 
C+ defined by: 



Imzl 



< 2, 1 Imzp > Scj^a^ X max(2dV^ax,da^ax)|- 



Lemma 6.6. Let T and T be given by Theorem 2.4. For every z €T>, 
there exists constants and (possibly depending on z) such that 



E 



(6.11) 



^T.-Riz)-^TrTiz] 



2+e/2 



< 



E 



-TTR{z)--Tif{z) 
n n 



2+E/2 



< 



and 



n 



l+£/4 ■ 
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Proof. We only prove the first inequality in (6.11). We use the identity 

R{z) - T{z) = R{z){T{z)-^ - Rizy'mz) 

= R{z){^{zy^ - B{zY^)T{z) - zR{z)A{^{z) - B{z))A^T{z). 
Therefore, j^TT{R{z) - T{z)) can be written as 



(6.12) 



1 Tr(i?(.) - T{z)) = 1 TTR{z){^{zr' - B-\z))T{z) 

-^TTR{z)A{i>{z)-B{z))A^T{z). 



The first term of the right-hand side of (6.12) is equal to 

l^f I 1 



lTri?(z)(^(z)-i - B-\z))T{z) = 1^ 



Let us first prove that the following control holds: 
(6.13) \{TR)ii{z)\< 



{TR)u{z). 



1 



I ImzP 

In fact, \{TR)ii{z)\ < \\Ti.{z)\\\\R.i\\ where Tj. (resp. R.i) denotes row num- 
ber i of T (resp. column number i of R). On the other hand, ||rj.(2;)|| < 
\im(z)\ ^ \im{z)\ Proposition 5.1 part 3. This proves (6.13). 

Minkowski's inequality leads to 



N 



TriRiz)i^izr'-B-'izmz)) 



2+e/2 



< 



1 



ImzP 



sup 



1 



1 



V'i(z) bi{z) 



2+e/2 



Similarly, the negative of the second term of the right-hand side of (6.12) 
can be written as 

n 

' Tr R{z)A{^iz) - B{z))A^T{z) = | Y.ii^,{z) - h,{z))a^T{z)R{z)^,. 



N 



Using again Minkowski's inequality and the fact that |ajr(z)i?(z)aj| < 



max 

|ImzP ' 



we get: 



-T,R{z)A{^{z)-B{z))A^T{z\ 



2+e/2 



< 



|z|(ia^ 



I ImzP 



-sup||V'j(2;) - &i(^;)||2+£/2- 
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As Tpjiz) — hj{z) is equal to — ^\z) ^^^'^^^ inequality 

hj{z)\\ < yields 

Uj{z)-b,{z)h+s/2< 
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hjiz) X 



1 



ipjiz) bj{z) 



2+S/2 



Gathering all the pieces together, we obtain 



(6.14) 



1tt{R{z)-T{z)) 



< 



1 



2+e/2 



Imzl 



■ sup 



1 



1 



+ 



\z\da:. 



2 

max 



ipi{z) hi{z) 
1 



I Imzl^ 



sup 



2+e/2 
1 



i)j{z) bj{z) 



2+e/2 



It is therefore sufficient to show that if z G then 



(6.15) 
(6.16) 



sup 

l<fc<Af 



sup 

i<e<n 



1 



1 



1 1 



tpi{z) be{z) 



2+e/2 



< 



2+e/2 



n 



where and do not depend on n and A^, but may depend on z. 

For this, we use (6.1) in the case where matrix U coincides with U = Di 
and constant Ki = cr^ax (recall that supj „ < c^ax)- Tr(-Di(5(2)) can 
be written as 



where supj ||ei(z) 



KS = K. 



l/(2+e/2). 



- TV DiQ{z) = - Tr DiR{z) + ei{z) 
n n 

||24-e/2 ^ ^) where A'f = T^^^^^"'"'^/^) (^e denote similarly 
. Using the very definitions of ipiiz) and of bi{z), we obtain 



1 



1 



ijjiiz) bi{z) n 



TiDi{R{z)-f{z)) + zei{z) 



Rewriting R{z) - f{z) as R{z) - f{z) = R{z){f{z)-^ - R{z)-^)f{z) and 
using similar arguments to what precedes [cf. (6.14)], we easily get 



1 



1 



(6.17) 



^i(^) bi{z) 



< 



Imzp 



2+e/2 



■ sup£ 



1 



'4)i{z) be{z) 



I kl ^max'^^max 

"T It I /I ^LIP 



|Imz|^ 



2+e/2 
1 



■ipkiz) bk{z) 



+ 



2+e/2 



n 
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Similarly, for each 1 < j < n, we have 
1 1 



iljjiz) bj{z) 



2+e/2 



(6.18) 



< 



max „ 

Tn-SUp 



1 



1 



_|_ 1^1 "max "'"•max 



1 



2+e/2 
1 



|Imz|4 "7 ^^(2) 5^(2) 

In order to simplify the notation, we put 
1 1 



+ 



2+e/2 



a 



n - sup 

l<A:<Af 



il^kiz) hk{z) 



2+e/2 



n - sup 

l<t<n 



1 



1 



2+e/2 



Equations (6.17) and (6.18) give immediately 



(6.19) 



\z\da 



I Im^l^ 

l^l'^'^ma; 

I ImzP 



I |2j 2 2 \-y\TC^ 
max- I ""max'^max ^, , I^K*-2 



I Imzp 



n 



max max - , 



Imzl^ 



n 



As zeV, we have '(imzl?'"'''' < ^- Therefore, 



2\z\do'i 



Imz 



2|z|K| 



n 



Plugging this inequality into (6.19), we obtain: 

('« on^ X ^ /^^I^P^^'-'^max , I^P'^'^max^max \ - , -^5 



Ima:!' 



Imzl' 



n 



for some constant > depending on z. It is then easy to check that 



|Imz|4 



+ ■ 



I |2j 2 2 



max max 



I Imzl 



< 1 



for z G P. Thus (6.20) implies that for z G P, q:„ < for some constant 
k^. Similarly, c^n < ^ and (6.15) and (6.16) are established. This, in turn, 
establishes (6.11). Proof of Lemma 6.6 is complete. □ 



6.3. Proof of Theorem 2.5. We are now in position to complete the proof 
of Theorem 2.5. We first remark that inequality (6.1) for U = I and inequal- 
ity (6.11) imply 

2+£/2 



(6.21) 



E 



lTrQ(.)-lTrr(.) 



< 



l+e/4 



n 
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for each z £T>. We consider a countable family {zk)kef>s with an accumula- 
tion point contained in a compact subset of P. Borel-Cantelli lemma and 
equation (6.21) imply that on a set Qi with probability one, -^Tr((5(zfc) — 
T{zk)) — > for each k. On Qi, the functions i— > Tr Q{z) and z ^ jj TvT{z) 
belong to 5(]R+). We denote by fniz) the function fniz) = jj-Tr{Q{z) — 
T{z)). On rJi, fn{z) is analytic on C — M+. Moreover, |/n(^)| < for each 
compact subset K (ZC — where 5k represents the distance between K 
and M^. By the normal family theorem (see, e.g., [21]), there exists a sub- 
sequence fuf. of /„ which converges uniformly on each compact subset of 
C - M+ to a function /* analytic on C - M+. But, f*{zk) = for each /c G N. 
This implies that /* is identically on C — M"*" , and that the entire sequence 
fn uniformly converges to on each compact subset of C — . Therefore, 
almost surely, Tr(Q(z) — T{z)) converges to for each z G C — M"*". Proof 
of Theorem 2.5 is complete. 

APPENDIX A: PROOF OF PROPOSITION 2.2 PART 3 

In this section, we prove Proposition 2.2 part 3. We first recall some results 
that concern the integral representations of some scalar complex functions 
analytic in the upper half plane = {z : Im(z) > 0}. 

A.l. The scalar case. These results can be found in Krein and Nudel- 
man's book [18] and therefore we adopt their notation in this section despite 
minor interferences with other notation in this article [beware in particular 
of the difference between class S functions below and Stieltjes transforms of 
probability measures denoted in the rest of the article by 5(M"^)]: 

Theorem A.l ([18], Theorem A.2). A function f{z) over C+ satisfies 
the following conditions: 

(i) f{z) is analytic over C'^ and 

(ii) Im/(z) >0 /or zgC+ 

if and only if it admits the (unique) representation 
(A.l) /(^)=a + 6z+ r f-^--^V(dO, 



where a G M, b>0 and fi is a positive measure such that 

fi{dt) 



l+t2 



< oo. 



The measure fi can be furthermore obtained by the Stieltjes inversion 
formula 

(A.2) \fi{{ti}) + \fi{{t2}) + K{ti,t2)) = - lim Im(/(t + ie)) dt. 
2 2 vr ejo Jfi 
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Different names are given to tfie functions f{z). In Krein and Nudelman's 
book [18], they are called class TZ functions. We will be particularly interested 
in the following subclass of these functions: A function is said to belong 
to class S if it belongs to class TZ and if it is furthermore analytic and 
nonnegative on the negative real axis (— cxd,0). 

Theorem A. 2 ([18], Theorem A.4). A function f{z) is in class S if and 
only if it admits the representation 

(A.3) ;(.) = e+r^, 

JO t — z 

where c > and v is a positive measure that satisfies 

uidt) 

< oo. 



l + t 

Class S functions can also be characterized by the following theorem: 

Theorem A.3 ([18], Theorem A. 5). A function f{z) is in class S if 
both f{z) and zf{z) are in class IZ. 

If f{z) belongs to 5, it also admits the representation 
(A,4) /(,) = „ + ,, + £(_L__^)^(,«). 

In the following, it is useful to recall the relationships between representa- 
tions (A.3) and (A.4). The intermediate steps of the proofs of [18], Theorem 
A.4 and [18], Theorem A. 5 give: 

Proposition A.4. The following relations hold: 

— ^ is carried by and fJ. = v, 

- b = 0, 

-a- /o+°° dii{t) >0 andc = a- /(+°° ^ dfi{t) > 0. 

We now address a partial generalization of these results to matrix- valued 
functions. 

A. 2. The matrix case. A matrix-valued function F{z) on is said to 
belong to class TZ if F{z) is analytic on C"^ and if ImF(z) > 0. Recall that 
matrix Im F{z) is defined as 

lmF{z) = ^^{F{z)-F*{z)). 

The generalization of Theorem A.l can be found, for example, in [8]: 
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Theorem A. 5 ([8], Theorem 5.4). An n x n function F{z) on C+ he- 
longs to class TZ if and only if it admits the representation 

(A.5) FW = .4 + ,B + £(^-^),W, 

where A is a Hermitian matrix, B >0 and fi is a matrix-valued nonnegative 
measure such that 

The proof is based on the corresponding result for the scalar case and the 
so-called polarization identity ([15], page 263) 

x*F{z)y = + yrF{z){x + y)-{x- yrF{z){x - y) 

+ i(x — \y)* F[z)[x — \y) — \{x + iy)* F{z){x — iy)). 

We are now in position to prove Proposition 2.2 part 3. It partly generalizes 
Theorem A. 2 to the matrix case: 



Theorem A. 6. If a matrix function F(z) satisfies: 
(i) F{z) and zF{z) are in class TZ 
then, it admits the representation 

fi{dt) 



(A.6) F{z)=C+r 

Jo 



t- z 



where C > and fx is a matrix-valued nonnegative measure carried by 
that satisfies 



Proof. Theorem A.6 could again be proved using the polarization iden- 
tity. We however provide the following shorter argument. Assume that F(z) 
and zF{z) are in class TZ. Then, F{z) can be written as (A.5). Let x S C". 
Then, x*F[z)x is scalar function which belongs to TZ. The representation 
given by (A.l) in Theorem A.l is given by 

roo / I t \ 
x*F(z)x = x*Ax + zx*Bx + f I ^ ]x* fix(dt). 

J~oo\t-Z l+t^J 

The quantity x*zF{z)x = zx*F{z)x also belongs to belongs to TZ. Therefore, 
it can be written as 

u^{dt) 



re 

x*F{z)x = c^+ / 
Jo 



t-z ' 
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where > and where Vx is a positive measure for which 

— < oo. 

1 + t 

Using Proposition A. 4, we get immediately that: 

- X* [XX is carried by M"*" and Vx = x* fix, 

- x*Bx = 0, 



x*Ax-J^j^ix*fix){dt)>0. 



The first item imphes that fj, is carried by and that Tv{J^ T^^) ^ 
The second item imphes that B = 0. As t{l + 1^)-! < 2(1 + ty^ for t > 0, 
the finiteness of Tt{J^{1 + 1)^^ ft{dt)) imphes 

Therefore F{z) can be written as 

F{z)=A- -^fj,{dt)+ 

Jo 1 + i Jo t - z 

The third item imphes that matrix C = A — /q°° j^/J'idt) is nonnegative. 
This completes the proof. □ 

We finally note that Theorem A. 6 can be found in several papers without 
proof. See, for example, [3], pages 64-65. 

APPENDIX B: PROOF OF PROPOSITION 6.4 

In the sequel, K will denote a bounding constant that does not depend 
on A^. Its value might change from one inequality to another. 

We first recall a useful result on quadratic forms associated with random 
matrices. 

A LEMMA FROM Bai AND SiLVERSTEiN (Lemma 2.7 in [1]). Let x = 
{Xi, . . . ,Xn)^ be a vector where the Xi 's are centered i.i.d. random variables 
with unit variance. Let C be an nx n deterministic complex matrix. Then, 
for any p>2, 

(B.l) E|x^Cx-T¥C|P<i^p((E|A:i|^TrC7C7*)P/2+E|A:ipPTr(C7C7*)P/2). 
B.l. Evaluation of Z^, and Z^. 



Proposition B.l. The random variables (5j, 7j and Pj satisfy 
sup E|/3,f +^ < sup E|7,f < J, sup E|/3,f < ^. 

where K is a constant that does not depend on n. 



(B.2) 
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Proof. Let us first establish the inequahty for (3j. We recaU that y^ can 
be written as yj = ■^D^-^'^-x.j where we recah that Xj = (-'^f j, • • • i^]v j)"^- 
Let V = Dy^Q^^^Aj. We can thus write (3j = ■^v'^Xj. We then have 

E[|aJg(^-)y.r+>] 

= ;^E[|xJvv-x/-/V] 
r,2+e/2 

< f__(E[|xJvv^x, - T¥ vv^|2+-/V] + I Tr vv^|2+-/2) 

+ E\Xu\^+' T¥(vv^(vv^)*)i+"/^ + II vf 

where (a) follows from the independence of v and Xj and from Lemma 2.7 in 
[1]. By noticing that ||v|| < ||L>]/^||sp||Q(j')||sp||aj|| < ^f^^ and by using 
Assumption A-1, we have the desired result. 

Let us now establish the inequality for jj. The random variable 7j can be 
written 

= yJi?(j)t/Q(^')aj + yJ{R - R^^^)UQ^^^a.j 

= yJi?(jVQ(^')aj + yji?(-'')(i?(-'')"' - R^^)RUQ^^^aij 

= yJi?(^-)C/g(^-)a, 

+ yJii(^')(S(j')"' - B-^ - - ABA^))RUQ'^^'>8Lj 

= 7i,i+7i,2-7i,3 + 7j,4, 

where 

7,-i = yJi2(jVQ(^)a,-, 
^^• 2 = yJi?(j)(B(j)"' - B~^)RUQ^i^a.j, 



(B.3) 



7i,3 = ^Y.^h^P - k)yjR^'^^i4RUQ^'^Sij, 
7^.4 = zbjyjR^^^aLjajRUQ'-^^aLj. 
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Beginning with the term jj^i, let v = R^^^UQ^^^aj. Recahing that yj = 
-^Dy^ X Xj and using the independence of v and y^, a standard calcu- 
lation leads to 

IE|7i,i|^<%^(EXi\+3)E||vf , 



where < (||i?(^')||sp||f/||sp||Q(^')||sp||a,||)^ < ^^f. This yields 
(B.4) E|7,,ir < ^^f^^iEXf, + 3)^^. 

Let us now consider 7^,2 = yJ-R^^'H^^^^"' - B-^)RUQ^^^aj . We have 

\%2\ < iii?(^')|i3pp(^')" -i?-lspii«iispnispiig(^')|ispiia,iiiiy,ii, 

(B.5) 

We now prove that 

(B.6) \\B^^^''-B-\^< 

By applying to (S^S — zl)~^ the matrix inversion lemma (see [16], Ap- 
pendix A, see also [15], Section 0.7.4), we obtain 





Z 


n 


Imz 



where 

^i,3 



lTrA^^')(S(^-)V^-)-z/)-\ 



n -y - 2^J(E0-)E0')^ - ziy^ij 



1 ^ 



n 



z — zt 

2 



-z-z^j(so-)sor -ziy^^j 



Definitions (2.7) and (6.2) yield 

-(T>(L)p^(s(j')' S(J') - ziy^) - TV(A(5]^S - z/)~^)). 



1 1 Z,_,~Q) u,T^ 



h^p{z) h{z) n 
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We have 

, 1 



Fi,2 = - 
n 



(B.7) 



n 



where (a) follows from a singular value decomposition of S^-'). In fact, let 
= ^^-^ f£U^vJ be a singular value decomposition of S'--'^ with i/^, U£, 
and V£ as singular value, left singular vector, and right singular vector re- 
spectively. Then, 



2| Tt |2 



1=1 1^^-^ 



Im(z|J(S(^)sO)^ - = Im(z) ^ ^ ' ^ 



1 



2 

which yields (B.7). Furthermore, one has |a;j.3| < ■ Thus, 



(B.8) 



< 



2cr2 


Z 


n 


Imz 



which yields (B.6). Plugging this into (B.5), we obtain 

n|Im# ll^^-ll- 
Finally, since E||yjf < a^^^d'^ (E\Xn\^ + 1), we have 

TC^S^U &'4„4 I |4 1 

(B.9) E|7,,/ < "T;^^Ct" (1 + lEl^nh oc ^ 

for N large enough. 
We now deal with 



7i,3 

£=1 



A rough control yields 



l7i,3| < |z|a^^||i2||,p||C/||,p||QW||,p^|6? -6,||yJ/?0-)a, 



^=1 
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,2 



We have 

hf{z)-Uz) 

= zb'f \z)be(z)-TT(DMJ^T.^ - zlV^ - (S^^'^S^^')'^ - z/)"^)). 
n 

Since 6^ and 6^"'' belong to 5(M~'") by Proposition 2.3 part 4(b), their absolute 
values are bounded above by | Imz|~^. By Lemma 2.6 in [25], we have 



(B.IO) |TrL>^((SS^ - ziy^ - (S(^')S(^')^ - < ^ 



2 

max 



I Imz| 

As a result, we obtain 

(B.ll) \hf{z)-b,{z)\<^ff^^. 
Thus, 



By Minkowski's inequality. 

The terms E|yji?(j)af |^ can be handled in the same way as for 7j^i. Thus, 
we obtain 

j4 I 18^12 12 7^4 -1 

(B.12) E|7,,3r < T2T ' (^l^nl' + 3) oc -. 

The last term 7j^4 = zbj-yj R^^^ ajaj RU Q^^^ aj satisfies 

Hence 

IE|7,,4r < 77:^7U2 aLxi^fE|yJi?(^')a 

(B.13) 




|4 



^Z.Ktai.Am^i\' + 3)z2 „2 



1 1 
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where (a) follows from the fact that the term K\yJ R^^^ajl^ can be han- 
dled as 7j-i. Gathering (B.4), (B.9), (B.12), (B.13) and using Minkowski's 
inequality, we obtain the desired inequality for 'jj. 

Using similar arguments, one can prove the same inequality for Pj. Propo- 
sition B.l is proved. □ 



The term Zf . Using (6.7), we have 



N ^^aj ~ Nllmzl ^ 
j=i J I I j=i 

Minkowski's inequality and Proposition B.l yield 



EN 



n|4 



< 



n\z\ 



N\ lmz\ 



K 



1 



2 2 ■ 



The term Z2 ■ We can write Z2 = jj Z)j=i ^j,2,i — Z]j=i ^j,2,2, where 
Zj,2,i = (1 - zhj{z)pj)(3j/aj and Z^- 2,2 = zbj{z)Pjpj/aj. 



We have \pj\ < \\Q 



-^)llsp||3-j|| < a4,ax/l I™-^!- As a consequence. 



|Z. 



'j,2,l 



< 



Imz 



(Imz) 



Proposition B.l yields: El^l"^ < K/n^. Therefore, E\Zj^2,i\^ < K/n'^ and 
Minkowski's inequality implies that ¥,\j^J2]'=i^j,2,i\'^ ^ K/v?. In particu- 
lar, E|i X;"=i ^i,2,iP+^/^ < K/n^+^'l'^ for e < 4 by Lyapunov's inequality. 
Let us consider Zj ^2,2- By the Cauchy-Schwarz inequality. 



mz. 



3,2,2 



2+e/2 ^ 



Imz 



4+£ 



We have 



|4+e 



< EC 



I 1 1 sp 1 1 1 1 sp 



E|4f+=)i/2(E|/3^f+-)i/2. 

Ily. 



'J I ; v^iMji 



(i)| 



< 



4+e 



(Imz)2 

By Minkowski's inequality, we have 

.4+e 



|4+£ 



E||yj 



1 4+e 



< 



AT 



2+e/2 



2+e/2 



n 



E ^X, 



2 



(B.14) 



^ m=l 



< 



CT- 



4+e 



n 



2+e/2 
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Therefore, E\(3j\'^+'^ < K. 

As E\(3j\^+^ < K/n^+^/^ by Proposition B.l, we have E| ^ Y.]j=i ^^,2,2^+"/^ < 
K/n^^^/^, which yields the desired result. 



The term Z^. We have 



1^3 



-. n ^ 



iV ^ a, 
j=i J 



„2 1^12 n 



N\lmz\3 ^ 

Therefore, Proposition B.l yields the desired result. 



Em- 



B.2. Evaluation of Z^- We rely here on Lemma 2.7 in [1] on the quadratic 
forms recalled at the beginning of this Appendix. Write = J2^=i ■2j-,4 
with 

Zj^A = (1 + zbj{z){l + ujj))pj/aj. 
We write ooj = TrDjQ{z) + ej^i + ej^2, where 

Ej^i =u;j-- Tr DjQ'-^^ (z) and e^- 2 = - Tr Dj (Q^^) (z) - Q{z)). 
n n 

Since ujj is equal to D]''^Q^^\z)D]''^^j, where = (Xf^, . . . , X]^ 
(B.l) yields 



< 



Ar^4 X l+£/4 /Vn-4+£ 
^n|4^'"max \ , Tt7lV" |4+e max 



TTji I 14-''" max i i TP I 1 4+e max 



n2+£/2 

Thanks to Assumption A-1, we therefore have 

By Lemma 2.6 in [25], we have 
(B.16) \ejo\ < 



2 

"^max 



n\ Imzl 
Therefore, 

\^jA = kPi^j(2;)/ajlki,i + ej,2| 
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where the last inequahty follows from \pj\ < \\Q^^\z)\\sp\\Sij\\^ < p^fff- Gath- 
ering (B.15) and (B.16), we obtain via Minkowski's inequality the desired 
result. 

B.3. Evaluation of Z^. Recall that ^5 = ^ E"=i ^ - f Tr(5^^ + zI)RUQ. 
We first prove that 

1 1 " 1 

(B.17) — Tv{B-^ + zI)RUQ = - - I] zqjj{z)—Tv DjRUQ. 

i=i 

This follows from the definition of bi{z). We have r-y-y + z = —zTvDiQ, thus 



B"^ + zl = diagf --Tr DiQ; l<i<N 



n 



- diag {'n);l<i<N 



n 

Z 



-E9jidiag(4(n);l<i<iV) 



n 

Z 



n . , 



which yields (B.17). Since ^ = —zqjj{z), we have 
1 n - 1 " 1 

,B.8) -V 

1 f 1 1 

We now study the asymptotic behavior of loj — jjTi{DjRUQ). Since 
= yjR{R^^^~' - i?-i)i?(j')[/Q(jVj 

= yjR{B^j'>~^ - zA'^^^B'^^^A^j^^ - B'^ + z^i^A^)i?(j')f/Q(jVj, 

we have 

^j-j^Tr{DjR{z)UQ{z))^X„ 
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where 

Xj = Xj,i + Xj,2 + Xj,3 + XjA 



with 



Xj,i = yjR^'^UQ^^^yj - ^ Tt{D,R{z)UQ{z)), 



Xj,2 



'J 

n 



Xj,3 = -zJ2{bf^ - 6,)yJi?a,aJ/2(^)c/Q(^)y„ 

XjA = zbjyjRa^SijR^^^UQ^^^yj. 

As usual we choose e > that satisfies Assumption A-1. We first handle 
Xj,i- Using the same arguments as those in Section B.2 to handle 

io, - - TtDjQ = -yjQ(^\z)y, - - Tr DjQ, 
n n n 

we obtain 

(B.19) E|xi,iP+'/^< ^ 



„i+e/4- 

The random variable Xj,2 satisfies 

\x,M < Pllsp||i?(^')" -i?-^|lspp(^''llspnispllQ(^')|lsplly.f , 

- |Imz|3" "^P"^^" - n|Imz|4 H^^H ' 

where (a) follows from (B.6). As a result, E|xj,2p^^^^ < j^2+e/2 E\\yj\ 
Using (B.14) we obtain 

(B.20) E|xj,2 ^ 



4+e 



Consider now the random variable Xj 3- Using the upper bound (B.ll) for 

(i) 



\bi — be\ , we have 



cr2 |z|2 " 



Minkowski and Cauchy-Schwarz inequalities yield 



(B.21) ||x,,3||2+e/2< 



" max I ^ I 

n\ Imzp 



X:(||yJi?a,||4+.||afi?(^-)C/Q(^Vill4+el 
.1=1 
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Consider first 

E|yJi?a,|4+-<E(p||^+-||a,r+-||y,r+-) 
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(B.22) 



.4+e 



< 



|4+e 



(a) 

<K, 



|Imz|4+e 

where (a) follows from (B.14). Let us now consider the term KlaJ R^^^UQ^^^ x 
yjl^^^ = ;^2T?72lE|v^xj|''+= with V = L>]/^Q(-')^C/i?(j)^a^ By a series of in- 
equalities similar to B.2, we obtain E|af i?(j)[/Q(j)yil^^^ < i^/n^+^Z^. In 
conclusion, we obtain 

(B.23) - - ^ 



Finally, the variable XjA can be handled in the same way 



(B.24) 



< 



< 



Imz 
K 



2+E/2 



(E|yJi2a,f+^)V2(E|aJi?(^VQ(^V,f 



„l+£/4- 

In order to finish the proof, it remains to gather equations (B.19), (B.20), 
(B.23) and (B.24) to get 

IE|Xi,i + Xj,2 + Xi,3 + XjA?-^''^ < 
Plugging this into (B.18) and using Minkowski's inequality, we obtain 

n 2+e/2 



E 



< 



K 



l+e/4 ' 



n 



which is the desired result. 



APPENDIX C: PROOF OF THEOREM 4.1 

C.l. Proof of the convergence of Cn{(T^) — Cn{cr^) to zero. The proof 
of the convergence relies on a dominated convergence argument. Denote by 
TTn the probability measure whose Stieltjes transform is -^TrT„(2;); sim- 
ilarly, denote by P„ the probability measure whose Stieltjes transform is 
^TVQ„(z): 



^TrT„(z) 



^ and 1tVQ„(z) 
X — z N 



X-z 



The following estimates will be useful. 
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Lemma C.l. The following equalities hold true 
(C.l) E 



(C.2) 



POO 1 1 

XTTnidX) 



l<i<N 
1<J<'^ 



Nn ^ N 

l<i<N 



In particular, 

(C.3) supffi r XFnidX)) = supf r X-KnidX)) < ^Lx + aLx- 



Proof. First notice that 



E / AP„(dA) = — ETVSS^, 
Jo 



which yields immediately equality (C.l). We now compute /q°° A7r„((iA). We 
first prove that 



(C.4) 
Compute 



;>oo 

/ XiTnidX) = lim Re 
Jo y-"^ 



iy^Trrn(i2/) + 1 



+ 1 



Hence, 



Re 



-iy(iy— Trr„(iy) + 1 



= iyiTn{dX) 
X-iy 

X^TTnjdX) 2 p A^„(dA) 



y 



2 X-KnjdX) 
A2 + y2 ■ 



The monotone convergence theorem yields limj^^oo 2/^ /o°° '^'^J'^^^^ = /o°° A7r„(dA). 
Equation (C.4) is proved. We now prove 

(C.5) lim -iy(iyr„(iy)+I^) = -diag(TrA,l<«<iV)+^^^. 

The mere definition of T„ yields 

T„(z) = {^-^ - zA^A^)-^ = {In - z^A^A^)'^^. 
Using the fact that 
{In - z^A^A^y'^ = In + z'^A^A'^ + {z'^A^A'^f{lN - z'^A^A^)''^, 
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we get that 

Tn{z) = * + z'^A^A^'^ + (z'^A^A^fTniz). 
We now compute 

-z{zTn{z) +In) 

(C.6) 

= -z(z^(z) + In) - z-^Az^A^z^ - z^Az-^A^z^Az-^A^Tniz). 

In order to compute the hmit of the previous expression when z = iy,y ^ oo, 
recah that hm— = lim— iyT„(i?/) =In and hm— iy^'(iy) =hm— iyx 
Tn{iy) = In whenever y — > cxd by (5.2) and (5.4). 

Let us first consider the first term on the right-hand side of (C.6). 

-iyiiy^iiy) + In) = diag( - - - iy 



■ diag 
1 



l + (l/n)T\-Ar„(iy) 

-{iy/n)TtDifn{iy) 
l + (l/n)TVAr„(iy) 

diag(Tr A, l<i<N). 



For the second term on the right-hand side of (C.6) we have 
{-iy^{iy))A{-iy^{iy))A^{-iy^{iy)) AA^ . 

The third term clearly converges to because Tn{iy) — > when y +oo. 
Equation (C.5) is established. This limit immediately yields 

Ji^m -i!,(i!,lTVT„(ii,) + l) = T^^IZ^-i + ^-^-^^ 

l<j<" 

Equating this equation with (C.4) implies that / A7r„((iA) is finite and gives 
its value. Therefore, equation (C.2) is proved. The inequality (C.3) follows 
immediately from (C.l) and (C.2). Proof of Lemma C.l is complete. □ 

We are now in position to prove that Cn(cr^) — C„((T^) 0. Recall that 

Cn{c7^) = ^ET¥Q„(-c^)) duj. 

The dominated convergence theorem together with Theorem 2.5 yield 

Vl^>0 --^^T,Tn{-Lo)-(--^¥.TTQn{-uo)\ ^0. 

OJ N \UJ N J n->oo 
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Moreover, 
1 



-UJ] 



< 



1 1 

CO /I 





Trr„(-(j) 
1 



+ 



1 1 



< 



X + io 



TTn{dX) 



+ 



E 



1 



1 



\uj X + (jj 



FnidX) 



+ 



E/o°°APn(dA) 



2 2 
^ 2 '^max ~^ ^mi 



which is integrable in lo over ((t2,oo). Therefore the dominated convergence 
theorem yields Cn{cr'^) — Cn{o''^) and the first part of Theorem 4.1 is 
proved. 

C.2. Proof of formula (4.3): some preparation. Performing the change 
of variable 7 = in (4.2) yields to the formula 



- 1 Tr -r 

7 V iV 7 



One can check that the integrand is continuous in zero. In fact, Lemma C.l 
yields 

^T^nidX) 



If 1 1/1 

- 1-— Tr-T 

TV N 7 V 7 



1 / 1 

7 V 7(^ + 7^"^) 



XTln{dX) 



XlTnidX) < 00, 



I + A7 7-^0, 

where 7r„ is the probability measure associated to the Stieltjes transform 
-^TrT(2:) We thus introduce slightly different notation than in the rest of 
the paper. These notations appear to be more convenient in the forthcoming 
computations. We denote by 



1/1 

5(7) = -t' 

7 



5(7) = -f 
7 



%(7) 







7 V 

















G = diag((9i,l <i < Af), 
e = diag(6ij,l< j<n). 



These notations yield 

5(7) = (9-^7) +7^e(7)^^)-\ 



%) = (6-^(7) + 7^^e(7)^)-^ 
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l + (7/n)TrA5(7)' 1 + (7/n) Tr Z?,-5(7) ' 

The general strategy to establish formula (4.3) is to write ^(1 — Tr 5(7)) as 
the derivative of some well-identified function of 7. The following quantities 
will also be of help: 

A = - diag(TYm,-5, 1 < i < n)), A = - diag(TY(A5, l<i<N)). 
n n 

In the sequel, we use both f'{'j) and ^(7) for the derivative of /. 

Lemma C.2. The following equality holds true: 
(C.7) SAQ = OAS. 

In particular, 

(C.8) Tr AQA'^S = Tr A^eAS, 

(C.9) Tr^G'^^5 = Tre'yl^e^5G-^ 

Proof. After elementary matrix manipulations (see also [15], Section 
0.7.4), we have 

(/ + uv)-^ = I-U{I + VU)~^V, 
which yields afterward 

(/ + uv)-^u = U{I + VU)'^. 

Let U = Q^I'^A and V = -leA^Q^/"^. Then 

(/ + 7eV2^eAT0i/2)-i0i/2^ ^ @y^A{i + -fSA'^QA)-^ 

^ SAQ = QAS, 

which is the desired result. Multiplying by , we obtain SAQ = A^QAS 
which yields (C.8). Multiplying to the left by Q' A!'^ and to the right by G"!, 
we obtain Q'A^SA = Q'A^eASQ-^ which yields (C.9). □ 

Lemma C.3. Denote by 

IM = -TttA'^Q'AS + — TtAQA^S'. 

N N 

Then the following equality holds true: 



(C.IO) 7(7) = ^ (^^ Tr(7A6^^ S) j - ^ (^^ log det 6(6-^ + 7^^ 6^) ) . 
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In order to prove Lemma C.3, we shall rely on a differentiation formula. 
Let A = A{x) he an N X N matrix, then the equality 

(C.ll) log det A(x) = Tr A-\x)^A(x) 

dx ax 

holds true in the case where A is Hermitian or A is the product of two 
Hermitian matrices, that is, A{x) = B{x)C{x). We provide here a short 
proof (in the case of a general square matrix, one can refer to [14], Section 
15). 

Consider first the case where A is Hermitian and write A = U AU* where 
A(x) = diag(Ai(x); 1 < i < iV) and UU* = [/*[/ = /. We have 

On the other hand, 

TvA~^-^A = TrUA^'^U*(U'AU* + UA'U* + UAU*') 
dx 

= Tr A-^A' + TrUA-^U*U'AU* + Tr UA-'^U*UAU*' 

= E #t + Tr(f/*f/' + U*'U) E 
~i Ai(x) ^ Xi{x) 

where (a) follows from the fact that U* U' + U*' U = which is obtained by 
differentiating U*U = I. 

In the case where A{x) = B{x)C{x) with B and C Hermitian, we have 

log <\etBix)C{x) = log det Bix) + log det Cix) 
dx dx dx 

= TiB-^B' + TvC-^C' 
= Ti B-^B'CC^^ + Ti C^^B-^BC 
= TtC~^B~^{B'C + BC') = TtA^^A', 
which is the expected result. We are now in position to prove Lemma C.3. 

Proof of Lemma C.3. We differentiate the expression 

= ^ Tr(AGA^5) + Tt{AQ'A^S) + ^ Tt{AQA^S') 
- ^ TV[(eA^GA)(/„ + jQA'^QA)-^] 



+ ^ Tr{Ae'A^S) - ^ Ti-[(e'^^e^)(I„ + 7GA^eA)-i] 
/(t) + ^ Tr{Ae'A^S) - ^ Tr[{Q'A^QA)iIn + jQA^^QA)''] 



DETERMINISTIC EQUIVALENTS FOR RANDOM MATRICES 53 

- ^ Tv[(e'^^eA)(/„ + 7e^^e^)-^] 

- ^ T¥[(eA^e'^)(/„ + jeA^QA)-'] 

/(7) + ^ Tr{AQA^S) - 1 TriA^QAS) 

^ Tr(Ae'A^S) - ^ 

N ^ ' N 

where (a) follows from the fact that {In + ^QA'^@A)~'^ = 56"\ (b) follows 
from (C.8) and (c) follows from (C.9). □ 

C.3. Proof of formula (4.3): the meiin system of equations. We are now 

in position to prove the following: 

Lemma C.4. Denote by 

J(7) = ^(l-^Tr5(7) 

Ji (7) = ^ Tr A'S + ^ Tr Ae'A'^S, 

^2(7) = ^ Tr A'S + ^ Tr A'^Q'AS. 
Then the following system of equations holds true: 

(C.12) A(^liogdet(e-i + 7AeA^)) = J(7) + ^i(7), 
(C.13) A(^^iogdet(e-i + 7A^e^)) = J(7) + J2(7), 



(C.14) 



A. (^^Tr(A + AQA'')S ) = 7(7) + Ji (7) 



+ ^Tr(A + ^e^^)5'. 



Proof. We first give equivalent formulations for the term ^(7): 
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^7 Tr((e-i - In)S + ^AQA^S) 



(C.15) = ^TiAS + ^TiAQA^S, 

(C.16) = ^Tr AS + ^Tr OAS, 

where (a) follows from (C.8) and from the fact that Tr AS* = Tr A5. Consider 
now 

log detre-^ + 7^6^^)^ = ^ Tr(e~^' + -fAQ'A^ + AeA^)S. 
d'y \N J N 

Easy computation yields Q^^' = A + 7A' and the previous equality becomes 

^ (^1 log det(e-i + 7^eA^)) = J(7) + Ji (7), 

where (C.15) has been used to identify J. Equation (C.12) is proved. One 
can prove similarly (C.13) by using (C.16). We now compute 

= ^ Tr A'S + ^TrAe'A^S +^TtAS + ^Tt AeA^S 

N N N N 

+ ^T]:{A + AQA^)S' 

J(7) + Ml) + +^Tr(A + A@A^)S', 
where (a) follows from (C.15). Equation (C.14) is proved. □ 

C.4. Proof of formula (4.3): end of the proof. Eliminating 7(7) between 
(C.13) and (C.14), we end up with 

±{1.Tt{A + A@A'')S^ 



Ji (7) + ^ log det(e-i + jA^eA) 
- ^ Tr A'S - T¥ A^e'AS + ^ Ti{A + AQA^)S' 



Since — ;^ Tr A'S + Tr A5' = 0, we obtain 



(C.17) 



= Ji(7) + (^logdet(e-i + 7A^GA)) +/(7). 
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First use (C.IO) which expresses 1(7) as a derivative, then extract Ji(7) 
from (C.17) and plug it into (C.12). After some simphfications, we obtain 

J(7) = log det (9^^ + jAQA')^ ' _ (^^ TV A^J + (^^ log det 
^F'(7). 

Since C{a^) = Jq^'^ 7(7)^7, a mere integration yields C{a'^) = F{a^^) — 
lim-y_»o -P'(7)- It remains to check that lim^_»o-^(7) = to obtain (4.3). In 
order to compute the limit of -P as 7 goes to zero, it is sufficient to check 
the following limits: 

lim S{'~f) = Ij\f, lim0(7)=/7v and lim0(7) = I„. 

7^0 7^0 7^0 

Let us prove the first limit, 

7 V 7/ 7JA + 7^ 7-»o 

where (a) follows from Proposition 5.1 part 2. In order to compute the 
limit involving Q and 0, it is sufficient to note that ^1(7) = ^ilji{—''y~^), 
to interpret 'i/'i as a Stieltjes transform (cf. Proposition 5.1 part 4) and to 
perform the same computation as for ^T(— 7"^). One can compute similarly 

the limit of 0. Theorem 4.1 is proved. 
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